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Abstract 



An overview of the theory and phenomenology of hadrons and QCD is provided from a Dyson- 
Schwinger equation viewpoint. Following a discussion of the definition and realisation of light-quark 
confinement, the nonperturbative nature of the running mass in QCD and inferences from the gap 
equation relating to the radius of convergence for expansions of observables in the current-quark 
mass are described. Some exact results for pseudoscalar mesons are also highlighted, with details 
relating to the Ua(1) problem, and calculated masses of the lightest J = 0, 1 states are discussed. 
Studies of nucleon properties are recapitulated upon and illustrated: through a comparison of the 
ln-weighted ratios of Pauli and Dirac form factors for the neutron and proton; and a perspective on 
the contribution of quark orbital angular momentum to the spin of a nucleon at rest. Comments 
\ on prospects for the future of the study of quarks in hadrons and nuclei round out the contribution. 

> 
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2 ! 1 Introduction 

q ; 

{sj . In trying to elucidate the role of quarks in hadrons and nuclei one steps immediately into the domain 
! of relativistic quantum field theory where within the key phenomena can only be understood via non- 
^ perturbative methods. Two prime examples are: confinement, the empirical fact that quarks have 
' not hitherto been detected in isolation; and dynamical chiral symmetry breaking, which is responsible, 
r~j ■ amongst many other things, for the large mass splitting between parity partners in the spectrum of 
^ ■ light-quark hadrons, even though the relevant current-quark masses are small. Neither of these phenom- 
. ?3 \ ena i s apparent in QCD's Lagrangian and yet they play a dominant role in determining the observable 
characteristics of real-world QCD. The physics of hadrons is ruled by such emergent phenomena. 



2 Confinement 

In connection with confinement it is worth emphasising at the outset that the potential between 
infinitely-heavy quarks measured in numerical simulations of quenched lattice-regularised QCD - the 
so-called static potential - is simply not relevant to the question of light-quark confinement. In fact, 
it is quite likely a basic feature of QCD that a quantum mechanical potential between light-quarks is 
impossible to speak of because particle creation and annihilation effects are essentially nonperturbative. 

A perspective on confinement was laid out in Ref. [1]. Expressed simply, confinement can be related 
to the analytic properties of QCD's Schwinger functions, which are often loosely called Euclidean-space 
Green functions. For example, it can be read from the reconstruction theorem that the only Schwinger 
functions which can be associated with expectation values in the Hilbert space of observables; namely, 
the set of measurable expectation values, are those that satisfy the axiom of reflection positivity [2] . This 
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is an extremely tight constraint. It can be shown to require as a necessary condition that the Fourier 
transform of the momentum-space Schwinger function is a positive-definite function of its arguments. 
However, that is not sufficient. 

In relation to 2-point Schwinger functions, which are those connected with the propagators of ele- 
mentary excitations in QCD, the axiom of reflection positivity is satisfied if, and only if, the Schwinger 
function possesses a Kallen-Lehmann representation. This statement is most easily illustrated for a 
scalar field, in which case it means that one can write the scalar-field's 2-point function in the form 1 

S( P 2 ) = f°cfc-£to P (<r) > V, > 0. (1) 
Jo V + T 

The spectral density for a non- interacting scalar field of mass m is p(q) = 5(q — m). It is plain that no 
function S(p 2 ) which falls-off faster than 1/p 2 at large spacelike momenta can be expressed in this way. 

An appreciation of the importance of the axiom associated with reflection positivity has led to the 
formulation of a confinement test [3]. With a momentum-space Schwinger function, S(p), in hand, the 
first step is to calculate 

A(r) = jd'xj -0^e^- s S( P ), (2) 

which gives the configuration-space Schwinger function in the rest frame. 2 One then examines the 
properties of A(t). 

If an asymptotic state of mass m is associated with this Schwinger function, then 

A(r) r =°°-^e- m ^ (3) 
v ; 2m w 

i.e., the Schwinger function is positive definite and the mass of the asymptotic, propagating state is 
given by 3 

— lim — In A(t) = ra . (4) 

t^oo (XT 

If, on the other hand, A(r) calculated from a particular Schwinger function is not positive definite, then 
the axiom of reflection positivity is violated and the associated elementary excitation does not appear 
in the Hilbert space of observables. Thus the appearance of a least one zero in A(T) is a sufficient 
condition for confinement. It is a very clear signal. 
An exemplar is provided by 

a function for which with any particular choice of e > there exists a p\ > such that Vp 2 > p\ one 
has \S(p) — 1/p 2 1 < e. One calculates from Eq. (5): 



A(t) = — e ^ T [cos fir — sin fir] . (6) 

This model was proposed in Ref. [5] and argued therein to describe a confined gluon because it doesn't 
possess a pole on the real-p 2 axis. The argument for confinement is supported by the fact that A(t) 



1 A Euclidean metric will be used throughout. In concrete terms that means: for Dirac matrices, {7^,7^} = 2<$ M „, 
7^ = 7^; and a ■ b = Y^t=i a *^»- ^ timelike vector, p^, has p 2 < 0. Naturally, no theory consistent with causality can 
produce a Schwinger function with a pole at spacclike p 2 . 

2 The arguments and conclusions that follow can readily be adapted to the case of models or theories without a mass 
gap; i.e., that possess a nonperturbatively massless excitation. 

3 This picture is readily generalised to the case of a Schwinger function that describes a channel with more than one 
asymptotic state; i.e., a ground state plus excitations [4]. 
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is not positive definite and, indeed, exhibits damped oscillations about zero. Hence the excitation 
associated with this Schwinger function cannot be connected with an expectation value in the Hilbert 
space of observables. 

The existence of a zero in A(r) and its disappearance have been advocated as a means by which 
deconfinement can be studied in QCD at nonzero temperature and density [6, 7]. Deconfinement and 
chiral symmetry restoration are coincident in all self-consistent studies of concrete models of QCD that 
exhibit both phenomena in vacuum. This result appears to follow from the crucial role played by the 
in-medium evolution of the dressed-quark self-energy in both transitions. Further studies are needed to 
fully elucidate the nature and the probable simultaneity of these transitions in QCD. 

It was observed and emphasised in Ref. [1] that a violation of reflection positivity is a sufficient but 
not necessary condition for confinement. This criterion is nevertheless a powerful discriminating tool, 
as can be illustrated through a concrete example, which also serves to introduce the Dyson- Schwinger 
equations (DSEs). The DSEs provide a nonperturbative approach in the continuum for calculating 
and modelling the properties of Schwinger functions. They have long been used to develop insight into 
gauge field theories at both weak and strong coupling [8]. 

In connection with confinement and dynamical chiral symmetry breaking (DCSB) the DSE for the 
dressed-quark propagator is of fundamental importance. In QCD this equation takes the form 

S( P y 1 = t 1 -p + m + J -^-tfD^ip- 9)y7^(9)rS(?,p), (7) 

where at this point it is sufficient to work with the unrenormalised equation. In Eq. (7), D^ v and r° 
are the dressed-gluon 2-point function (propagator) and dressed-quark-gluon 3-point function (vertex); 
and m is a current-quark bare mass. Finally, the subject of the equation, S(p), is the dressed-quark 
propagator in vacuum. It can be written in a number of equivalent forms: 

S(p) = -iTpaytf) + *s(p 2 ) = , 7 . M(p2 1 ) + g(p2) = ry ■ p + M(p 2 ) ' (8) 

each of which involves just two scalar functions. In the last representation they are chosen to be the 
wave function renormalisation, Z(p 2 ), and the running mass function, M(p 2 ). 
Reference [9] introduced a simple Ansatz for the kernel of Eq. (7); namely, 



g 2 D llu (k) = (2n) i g6\k) 



rS(?,p) = y7„, (9) 



where Q is a mass 2 -scale. The ladder quark-antiquark interaction represented by Eqs. (9) is a constant 
in configuration space. Therefore the usual proof of the cluster decomposition property fails and hence 
this Ansatz defines a confining model. Furthermore, with m — 0, Eqs. (7) and (9) yield the following 
solution: 



A(p) ~ ||(i + > /r+W); v 2 >Q ' B(p) ~ \o ; p*>g ■ (10) 

When Eq. (2) is evaluated with as(p 2 ) constructed from Eqs. (10), one obtains 



AM«^f; (ii) 

i.e., a function with repeated zeros. Thus the interaction's violation of the cluster decomposition 
property is expressed in a fermion propagator that fails to satisfy the axiom of reflection positivity. 
Hence it describes a confined elementary excitation. 
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The story continues in a fascinating manner when one dresses the quark-gluon vertex. Studies 
attempting to elucidate the nature of p) are now underway using both continuum [10, 11, 12, 13, 14] 
and lattice-QCD techniques [15]. However, this is a relatively recent development. It was and remains 
common to employ an Ansatz for this 3-point Schwinger function that is constructed according to 
lessons learnt by studying Ward-Takahashi identities in gauge field theories [16, 17, 18]. One often used 
form is based on the vertex introduced in Ref. [19]; viz., 4 r^(£i,£ 2 ) = 4p T^i, £2), with 

iT^iM = t^A(£i,el)i, + (ei + e 2 ),[tli-(£i + e2)A A (£lf 2 ) + A B (£l£ 2 2 )], (12) 
Vf(%,® = \ mi) + f{£1)\ , ^f{£UD = ni \ I *P , (13) 

where F = A,B; viz., the scalar functions in Eq. (8). If one solves the gap equation, Eq. (7), using this 
vertex model and g 2 D flu from Eq. (9) one finds in the chiral limit 



< i\ p 2 -2g[l-exp(V/[2g])] ,2x1 2 
<W ) = ~ 4 , °s{P ) = T^jg exp(-p /[2Q\) . (14) 

It is apparent that dressing the vertex has markedly altered the form of the dressed-quark propagator; 
namely, ays, which each had a branch point with the bare vertex, are transformed by the dressed vertex 
into entire functions. Moreover, when Eq. (2) is evaluated with o~s(p 2 ) constructed from Eqs. (14), one 
obtains 

A(t) oc exp(-^r 2 /2) . (15) 
In contrast to Eq. (11), this is positive definite. However, it yields an effective mass via Eq. (4): 

m = lim m(r) = lim Qr = oo , (16) 

and infinitely massive states cannot propagate. 

These observations notwithstanding, confinement of the quark represented by the propagator in 
Eq. (14) is still signalled by a violation of the reflection positivity axiom. It is clear from the momentum 
space form of the propagator that the axiom is breached because both o~v t s fall faster than 1/p 2 at large 
spacelike p 2 . However, it can also be seen in another way, which is indispensable if the momentum-space 
form of the Schwinger function is not known. One can evaluate Eq. (2) with Eq. (1) to obtain 

A(r) = /°°*4r e ""- (17) 

Reflection positivity is satisfied by a configuration-space Schwinger function if, and only if, it can be 
expressed in this form, with p(q) as described in Eq. (1). Now the theorem of Ref. [21] can be employed 
to show that the function in Eq. (15) cannot be represented in this way and therefore contravenes the 
axiom. It is a fact that no function which falls faster than exp(— ar) at large r, where a is some constant, 
can be expressed through Eq. (17). 5 Amongst other things, it follows that an excitation described by 
an effective mass, m(r), which grows with Euclidean separation is necessarily absent from the spectrum 
of physical states because its Schwinger function violates the axiom of reflection positivity. This is a 
novel perspective. 



4 Improvements can be considered. However, modifications that vanish at t\ = £2 do not change the results which 
follow. The influence of more significant changes; for example, introducing a dependence on functions other than those 
which appear in the quark propagator, must be explored on a case-by-case basis. It is a simple truth that all we have at 
present are Ansdtze for the vertex and no modification will alter a simple, long established fact, soon to become apparent: 
the vertex can have a material impact on the expression of confinement and DCSB in QCD. 

5 Correspondence with W. N. Polyzou was instrumental in developing this aspect of the discussion. 
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With the information provided above one can identify a confined excitation; and hence, in the context 
of QCD, these tools have been employed to analyse propagators obtained both as solutions of truncated 
systems of DSEs and through numerical simulations of lattice- regularised QCD. Extant analyses, e.g., 
Ref. [22], suggest strongly that the gluon 2-point function in QCD is not reflection positive. The same 
has been argued for the dressed-quark propagator in QCD, e.g., Ref. [11]. However, as noted above, 
the detailed structure of the dressed-quark-gluon vertex is unknown and it can influence the analytic 
properties of S(p). Hence a firm conclusion on the dressed-quark propagator has not yet been reached 
(see, e.g., Ref. [23]). 

In relation to the studies just mentioned, it is noteworthy that any 2-point Schwinger function with 
an inflexion point at p 2 > must breach the axiom of reflection positivity. This can be seen directly 
from Eq. (1), which constrains the first derivative of reflection positive functions to be negative definite 
and the second derivative to be positive definite. Hence, it is often unnecessary to actually calculate 
the Fourier transform described in Eq. (2). A violation of reflection positivity can be determined by 
inspection of the pointwise behaviour of the momentum space Schwinger function. 

It should now be clear that we have a means of identifying whether or not a given Schwinger function 
can be associated with an expectation value in the Hilbert space of observables. In principle, one could 
calculate all QCD's Schwinger functions and demonstrate confinement of colour therewith. Plainly, 
however, that will be difficult in practice and nigh on impossible by brute force. A proof of confinement 
at some higher organisational level is necessary, and sought by many. 

There is a surviving hope that the study of gauge theories in less than four dimensions might be 
useful in this regard. Analogous to quenched QCD, quenched QED in three dimensions (two spacial, one 
temporal - QED3) is confining because it has a nonzero string tension [24]. The effect of unquenching; 
viz., allowing light fermions to influence the theory's dynamics, remains an important but unresolved 
question. It is conceivable that through its answer we might inform ourselves about aspects of the 
confinement problem in QCD. The early history of QED 3 studies is reviewed in Ref. [8] and their 
current status can be traced from Ref. [25]. 

In closing this section it is important to reiterate that confinement is a conjecture. In order to 
establish confinement as more than a contemporary empirical fact, it will be necessary to map out 
the nonperturbative structure of QCD's /^-function. 6 In order to accomplish that, one will need to 
determine the nonperturbative evolution of the renormalisation constants. That is tantamount to 
establishing rigorously that the theory exists. This fact emphasises the magnitude of the problem. 

Absent this, progress in hadron physics can still be made with well-constrained models that express 
key consequences of confinement. The feedback that this enables between experiment and theory further 
illuminates the essence of the confinement conjecture in relation to light-quarks. This can begin with 
the development and employment of a Poincare covariant framework that ensures the colour-singlet 
S'-matrix elements associated with physical processes do not possess production thresholds for light- 
quarks, in particular, and coloured elementary excitations in general. The notions described above 
have been immensely valuable in this regard. For example, they underly the DSE model introduced in 
Ref. [27], which has served as an archetype for numerous efficacious applications, as can be seen, e.g., 
from Refs. [28, 29, 30, 31]. 



6 Thc requirement that renormalisation preserve gauge invariance; viz., that the rcnormalised QCD Lagrangian be 
invariant under Bccchi-Rouct-Stora transformations, leads to the Slavnov- Taylor identities, which entail that all Schwinger 
functions which can be used to define the running coupling must yield the same result for that coupling. It follows that 
QCD possesses a unique /^-function. For a detailed discussion, see, e.g., Ref. [26]. 
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Figure 1: Dressed-quark mass function, M(p): solid curves - DSE results, obtained as 
explained in Refs. [11, 33], "data" - numerical simulations of unquenched lattice-QCD [34]. 
In this figure, adapted from Ref. [35], one observes the current-quark of perturbative QCD 
evolving into a constituent-quark as its momentum becomes smaller. The constituent-quark 
mass arises from a cloud of low-momentum gluons attaching themselves to the current- 
quark. This is dynamical chiral symmetry breaking: an essentially nonperturbative effect 
that generates a quark mass from nothing; namely, it occurs even in the chiral limit. 



3 Dynamical Chiral Symmetry Breaking 

Many statements of fact can be made in connection with this emergent phenomenon. For example, 
DCSB explains the origin of constituent-quark masses and underlies the success of chiral effective field 
theory. Understanding DCSB within QCD proceeds from the renormalised gap equation [32]: 

S^)- 1 = Z a (ry • p + m bm ) + Z x jf Vv(p - q)^l,S{q)Tl{q, P ), (18) 

where represents a Poincare invariant regularisation of the integral, with A the regularisation mass- 
scale, D^ v is the renormalised dressed-gluon propagator, T u is the renormalised dressed-quark-gluon 
vertex, and m bm is the quark's A-dependent bare current-mass. The vertex and quark wave-function 
renormalisation constants, Z lj2 (C 2 , A 2 ), depend on the gauge parameter. 

The form of the solution to Eq. (18) is the same as that written in Eq. (8) except that in QCD one 
must account for the renormalisation point dependence: 

S(p) = —27 ■ pa v (p\ C 2 ) + os(p\ C 2 ) = %i . M(p2j ^ - C2) = %1 Z p (P +mU ■ (19) 

It is important that the mass function, M(p 2 ) = B(p 2 , ( 2 )/A(p 2 , ( 2 ), illustrated in Fig. 1, is independent 
of the renormalisation point, (. The dressed propagator is obtained from Eq. (18) augmented by the 
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renormalisation condition 

S(p)- l \ p2=( . 2 =tl-p + m(( 2 ), (20) 

where m(( 2 ) is the running mass: 

Z 2 (( 2 , A 2 ) m bm (A) = Z 4 (C 2 , A 2 ) m(C 2 ) , (21) 

with the Lagrangian-mass renormalisation constant. In QCD the chiral limit is strictly defined by 
[32]: 

Z 2 (C 2 ,A 2 )m bm (A)^0,VA 2 >C 2 , (22) 

which states that the renormalisation-point-invariant current-quark mass m — 0. 

In perturbation theory it is impossible in the chiral limit to obtain M{p 2 ) 7^ 0: the generation of mass 
from nothing is an essentially nonperturbative phenomenon. On the other hand, it is a longstanding 
prediction of nonperturbative DSE studies that DCSB will occur so long as the integrated infrared 
strength possessed by the gap equation's kernel exceeds some critical value [8]. There are strong 
indications that this condition is satisfied in QCD [33, 34] (see Fig. 1). It follows that the quark-parton 
of QCD acquires a momentum-dependent mass function, which at infrared momenta is ~ 100-times 
larger than the current-quark mass. This effect owes primarily to a dense cloud of gluons that clothes 
a low-momentum quark [35]. It means that the Higgs mechanism is largely irrelevant to the bulk of 
normal matter in the universe. Instead the single most important mass generating mechanism for light- 
quark hadrons is the strong interaction effect of DCSB; e.g., one can identify it as being responsible for 
98% of a proton's mass. 

The vacuum quark condensate is the order parameter most commonly cited in connection with 
DCSB. For example, it occurs as a parameter at leading order in chiral perturbation theory. It is a 
basic fact that only in the chiral limit is it possible to unambiguously define the gauge invariant vacuum 
quark condensate in terms of S(p) [32, 36, 37]. That such a definition is possible at all emphasises that 
gauge covariant quantities contain gauge invariant information. Notwithstanding these truths, M(p 2 ) 
is a more fundamental tracer for DCSB. The condensate is only a small part of the information that 
the mass function contains, amongst other pieces are pointwise consequences of Goldstone's theorem 
[32], definitions and an explanation of the constituent-quark mass [35, 38, 39], and a related a-term 
[40]. 8 Harking back to Sect. 2, it is interesting to note that considered as a function of current-quark 
mass the dressed-quark propagator possesses a spectral representation [36], confinement and DCSB 
notwithstanding. 

That a nonzero quark mass function is generated in massless QCD entails that chiral symmetry is 
realised in the Nambu-Goldstone mode and the vacuum quark condensate is nonzero. The magnitude 
of the mass function in the infrared and order parameters contingent thereupon underly the validity of 
chiral effective theory as a tool for correlating low-energy observables in QCD. The question of whether 
M{p 2 ) has an expansion in current-quark mass around its chiral-limit value bears upon the radius of 
convergence for that perturbation theory. This question asks whether it is possible to write 

00 

M{p 2 ;m) = M{p 2 ;m = 0) + ^a n m n (23) 

n=l 

on a measurable domain of current-quark mass. It was found [37] that such an expansion exists for 
m < m rc and is absolutely convergent. The value of m rc can be reported as follows: 9 for a pion-like 

7 Owing to asymptotic freedom it is natural to fix the dressed-propagator to be free-particle- like at some large spacelike 
momentum scale, £. 

8 The current-quark mass dependence of both this cr-term and the massive-quark condensate defined unambiguously 
in Ref. [37] demonstrate that the essentially dynamical component of chiral symmetry breaking decreases with increasing 
current-quark mass. This is an important observation that is basic, e.g., to the validity of heavy-quark effective theory. 

9 NB. Irrespective of the current-mass of the other constituent, a pseudoscalar meson containing one current-quark 
whose mass exceeds m cr is never within the domain of absolute convergence. 
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meson constituted from a quark, /, with mass m rc and an equal-mass but different flavor antiquark, 
g, m Q gf = 0.45 GeV. Since physical observables, such as the leptonic decay constant, are expressed via 
M(p 2 ), it follows that a chiral expansion is meaningful only for (m^) 2 < 0.2 GeV 2 . Hence, it is only 
valid to employ chiral perturbation theory to fit and extrapolate results from numerical simulations of 
lattice-regularized QCD when the simulation parameters provide for m 2 < 0.2 GeV 2 . Lattice results at 
larger pion masses are not within the domain of convergence of chiral perturbation theory. 

4 Axial-vector Ward-Takahashi identity and mesons 

Dynamical chiral symmetry breaking is a fact in QCD and this amplifies the importance of the axial- 
vector Ward-Takahashi identity. The existence of a sensible truncation of the DSEs [41, 42] has enabled 
proof via that identity of a body of exact results for pseudoscalar mesons. They relate even to radial 
excitations and/or hybrids [43, 44, 45], and heavy-light [46, 47] and heavy-heavy mesons [39]. The 
results have been illustrated using a renormalisation-group-improved rainbow-ladder truncation [38, 48], 
which also provided a prediction of the electromagnetic pion form factor [49]. In addition, algebraic 
parametrisations of the dressed-quark propagators and meson bound-state amplitudes obtained from 
such studies continue to be useful. One example is a recent application to S-meson — > light-meson 
transition form-factors [50], the results of which can be useful in the analysis and correlation of the 
large body of data being accumulated at extant heavy-quark facilities, and thereby in probing the 
Standard Model. 

4.1 Neutral pseudoscalars and the rj' 

Implications for neutral pseudoscalar mesons have been elucidated [51]. In the general case the axial- 
vector Ward-Takahashi identity is written 

P^(k; P) = S-\k + )i lb F a + i lb T a S-\k_) - 2iM ab T b 5 (k; P) - A a (k; P) , (24) 

wherein: {jF a |a = 0, . . . , N 2 — 1} are the generators of U(Nj); the dressed-quark propagator S = 
diagfS^, Sd, S s , S c , Sb, . . .] is matrix- valued; -M(C) is the matrix of renormalised (running) current-quark 
masses and M ab = tr F [{J ra ,M}j ub ] , where the trace is over flavour indices. The inhomogeneous 
axial- vector vertex in Eq. (24) satisfies a Bethe-Salpeter equation; viz., 

[aft P)] tu = Z 2 [757^L + /"V(<?+)Il,(g; P)S(«Z-)U*£(?, k; P) , (25) 

Jq 

where: P is the total momentum and k, q are relative momenta, k± = k±P/2; r,. . . ,u represent colour, 
Dirac and flavour indices; and the quantity K(q, k; P) is the Bethe-Salpeter kernel. The pseudoscalar 
vertex satisfies an analogous equation driven by the inhomogeneity Z^T a 1 . 
The final term in Eq. (24) expresses the non-Abelian axial anomaly. It involves 

Au(k-P) = Jd 4 xd 4 ye l(k+ - x - k - y) N f (F°q(x) Q(0)q(y)), (26) 

wherein the matrix element represents an operator expectation value in full QCD and 

= ^ e^F^F^x) = d,K,{x) (27) 

is the topological charge density operator, with F° the gluon field strength tensor. It is fundamentally 
important that while Q(x) is gauge invariant, the associated Chern-Simons current, K^, is not. Thus 
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in QCD no physical state can couple to K^. Hence, physical states cannot provide a resolution of the 
so-called £/4(l)-P r oblem; i.e., they cannot play any role in ensuring that the rj' is not a Goldstone mode. 

It is plain and important that only A a=0 can be nonzero. Moreover, in considering the ?7A(l)-P r oblem 
one need only focus on the case A ^ because if that is false, then following Ref. [32] it is clear that 
the rj' is certainly a Goldstone mode. ^4° is a pseudoscalar vertex and can therefore be expressed 

A\k; P) = ^° 75 [i£ A (k; P) + 7 • PT A {k- P) +1 .kk- PQ A (k; P) + a^P u H A (k; P)} . (28) 

Equation (24) can now be used to derive a collection of chiral-limit, pointwise Goldberger-Treiman 
relations, important amongst which is the identity 

2f r , l E T]l (k; 0) = 2B (k 2 ) - £ A (k; 0) , (29) 

where B (k 2 ) is obtained in solving the chiral-limit gap equation. 
It is now plain that if 

£ A (k;0) = 2B (k 2 ), (30) 

then frfE^ik; 0) = 0. This being true, then the homogeneous Bethe-Salpeter equation for the rf does 
not possess a massless solution in the chiral limit. The converse is also true; namely, the absence of 
such a solution requires Eq. (30). Hence, Eq. (30) is a necessary and sufficient condition for the absence 
of a massless rj' bound-state. It is the chiral limit that is being discussed, in which case B (k 2 ) ^ if, 
and only if, chiral symmetry is dynamically broken. Hence, the absence of a massless rj' bound-state is 
only assured through the existence of an intimate connection between DCSB and an expectation value 
of the topological charge density. A relationship between the mechanism underlying DCSB and the 
absence of a ninth Goldstone boson was also discussed in Ref. [52]. 

Reference [51] derives a range of corollaries, amongst which are mass formulae for neutral pseu- 
doscalar mesons. These led, e.g., to a prediction of the manner by which the r( is split from the octet 
pseudoscalars by an amount that depends on QCD's topological susceptibility. That is most easily 
illustrated by considering the U(Nf) limit, in which all current-quark masses assume the single value 
m((). In this case one finds 

m^=rv + 2m(C)pJ(C), (31) 

where m v i is the meson's mass, 

/° P M = Z 2 tr / ^° 757m xrf (q\P), ip'iO = Z A tr [ ^° 75 Xr/ (q; P) , (32) 

J q J q 

with Xiri(k] P) = S(k + )r 7r .(k] P)S(k_) } and the residue of the v[ bound-state pole in ^4° of Eq. (26) is 

*V = ^iV> "* = <0|e|i/>. (33) 

Plainly, the rj' is split from the Goldstone modes so long as n v r ^ 0. 10 
It is argued [54, 55] that in QCD 

n v' ~ Ty^f , ( 34 ) 
and it can be seen to follow from the gap equation, the homogeneous BSE and Eqs. (32) that 

fi~VN e ~P$(Q. (35) 



10 A nonzero value of the topological susceptibility can be achieved through the coupling of a massless axial-vector 
gauge-field ghost to the Chern-Simons current, which does not appear in the particle spectrum of QCD because the 
current is not gauge invariant. This is a variant of the Kogut-Susskind mechanism [53]. 
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Table 1: Masses (GeV) of the lightest J = 0, 1 states produced by the rainbow-ladder DSE 
truncation of Refs. [38, 48] with the parameter values: uj = 0.4 GeV, uD = (0.72 GeV) 3 ; and 
current-quark masses m M) d(lGeV) = 5.45 MeV, m s (l GeV) = 125 MeV. The rainbow-ladder 
kernel gives ideal flavour mixing for all states, and hence the pure ss 0~ + state is an artefact 
which does not correspond to a physical state. A simple kernel that implements flavour 
mixing is described in Ref. [51]. The text explains in detail the results in this table, which 
is adapted from Ref. [62]. 
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Consider now Eq. (31) in the form 

^, = ^ + 2m(C)^P. (36) 

J v > J v ' 

The first term is zero in the limit N c — > oo while the second remains finite. Subsequently taking the 
chiral limit, m v > vanishes in the manner characteristic of all Goldstone modes. 11 These results are 
realised in the effective Lagrangian of Ref. [57] in a fashion that is consistent with all the constraints of 
the anomalous Ward identity. This is not true of the so-called 'tHooft determinant [58, 59, 60]. 

Reference [51] also presents an Ansatz for the Bethe-Salpeter kernel that enables illustration of the 
implications of the mass formulae, in particular mixing between neutral pseudoscalars. 12 Despite its 
simplicity, the model is elucidative and phenomenologically efficacious; e.g., it predicts r\-r{ mixing 
angles of ~ —15°, n -?] angles of ~ 1°, and a strong neutron-proton mass difference of 0.75 (m d — rh u ). 



4.2 Meson spectroscopy 

The use of DSEs to study meson phenomena is empowered by the existence noted above of a systematic, 
nonperturbative and symmetry-preserving truncation scheme [41, 42]. It means that exact results, 
such as those outlined heretofore, can be illustrated and predictions for experiment made with readily 
quantifiable errors. The renormalisation-group-improved rainbow-ladder truncation of the gap and 
Bethe-Salpeter equations introduced in Refs. [38, 48] has been widely employed in this endeavour. To 
exemplify that, in Table 1 calculated results are reported for the masses of the lightest J = 0, 1 states 
[62]. It is true in general that the truncation is accurate for the _+ and 1 light-quark meson ground 
states. In these channels it can be seen algebraically that contributions beyond rainbow-ladder largely 
cancel between themselves owing to Eq. (24) [12, 14, 42, 63]. The remaining columns in the table deserve 
special attention because they show clearly the path toward improvement. 

Terms beyond the rainbow-ladder truncation are known to add constructively in the ++ channel 
[64]. Hence the leading order truncation is a priori not expected to provide a good approximation. 
Further understanding is provided by an exploration of the contribution from two-pion intermediate 
states to the mass and width of this lowest-mass scalar. A rudimentary analysis shows that a realistic 
description is attainable therewith [40]; viz., it gives a pole position y/s a = (0.578 — z 0.311) GeV. This 
is not the end of the scalar story but it is a sensible path to follow, in particular because a QCD-level 
mechanism is precisely specified. Something analogous should be considered in connection with the ss 

n NB. One must take the limit N c — > oo before the chiral limit because the procedures do not commute [56]. 
12 Interplay between the non-Abelian anomaly and flavor symmetry breaking is also explored in Ref. [61]. 
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scalar state listed in Table 1. It is notable that in rainbow-ladder truncation it has been shown that at 
least up to the c-quark mass the ordering of meson masses is 0~ + (lS') < ++ (lS') < 0~ + (2S) < ++ (2S), 
where n — 1,2 denotes the ground state and first radial excitation, respectively [65]. 

Compared with experiment, the masses of the axial- vector mesons l +± are poorly described by the 
rainbow-ladder truncation: ~ 400 MeV of repulsion is missing from the kernel. A cruder model does 
better [66, 67]. The latter studies and a more recent analysis [68] indicate that at least part of the 
defect owes to the absence of spin-flip contributions at leading-order. Such contributions appear at all 
higher orders and are enhanced by the strongly dressed quark mass function. This is one of the ways 
that the meson spectrum can be used to probe the long-range part of the light-quark interaction and 
thereby to chart the nonperturbative behavior of QCD's /3-function. 

The last three columns describe systems with so-called exotic quantum numbers. Of course, these 
states are exotic only in the context of the naive constituent quark model. In QCD they correspond 
simply to interpolating fields with some gluon content and are easily accessible via the BSE [69]. 
Nonetheless, while the rainbow-ladder truncation binds in these channels, the shortcomings encountered 
in the 1 + channels are also evident here, for many of the same reasons. Reliable predictions for the 
masses of such states will only be obtained once improved kernels are developed. At the very least, 
one must have dependable predictions for axial-vector masses before drawing any conclusions about 
the so-called exotics. Since all these states lie within a domain on which DCSB is very relevant, a 
framework with a veracious description of that emergent phenomenon is essential. 

One might pose the question of whether, in the context of bound-state studies in which model 
assumptions are made regarding the nature of the long-range interaction between light quarks, anything 
is gained by working solely with Schwinger functions. This means, in part, constraining oneself to work 
only with information obtained from the DSEs at spacelike momenta. 13 According to a recent study 
[4] the answer is no. It analysed the capacity of Schwinger functions to yield information about bound 
states, and established that for the ground state in a given channel the mass and residue are accessible 
via rudimentary methods. However, simple methods cannot provide dependable information about 
more massive states in a given channel. Indeed, there is no easy way to extract such information. 
An approach based on a correlator matrix can be successful but only if the operators are carefully 
constructed so as to have large overlap with states of interest in a given channel, and statistical and 
systematic errors can be made small; viz., ~ 1%. While it is possible in principle to satisfy these 
constraints, doing so is labor intensive and time consuming. That is only justified in the absence of 
model-dependence. 

5 Nucleons 

Partly because attention to this sector has only recently increased, the current level of expertise in 
studying baryons is roughly the same as it was with mesons more than ten years ago; viz., model 
building and phenomenology. We are a little ahead of that game, however, because much has been 
learnt in meson applications; e.g., as indicated above, a veracious understanding of the structure of 
dressed- quarks and -gluons has been acquired. 

5.1 Faddeev equation 

In quantum field theory a nucleon appears as a pole in a six-point quark Green function. The pole's 
residue is proportional to the nucleon's Faddeev amplitude, which is obtained from a Poincare covari- 

13 Lattice-regularised QCD provides a background to this question. That approach is grounded on the Euclidean space 
functional integral. Schwinger functions are all that it can directly provide. Hence, lattice-regularised QCD can only be 
useful if methods are found so that the question can be answered in the affirmative. 
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ant Faddeev equation that adds-up all possible quantum field theoretical exchanges and interactions 
that can take place between three dressed-quarks. This is important because modern, high-luminosity 
experimental facilities employ large momentum transfer reactions; viz., Q 2 > Mfj where M N is the nu- 
cleon's mass, and hence a veracious understanding of contemporary data requires a Poincare covariant 
description of the nucleon. 

A tractable truncation of the Faddeev equation is based [70] on the observation that an interaction 
which describes mesons also generates diquark correlations in the colour-3 channel [71]. The dominant 
correlations for ground state octet and decuplet baryons are + and 1 + diquarks because, e.g.: the 
associated mass-scales are smaller than the baryons' masses [66, 72], namely (in GeV) - 

™>[ud\ + = 0.7 - 0.8 , m (uu ) i+ = m {ud)i+ = m {dd)i+ = 0.9 - 1.0 ; (37) 

and the electromagnetic size of these correlations, while larger than the pion, is less than that of the 
proton [73] - 

f[ud\ 0+ ~0.7fm, r (ud)i+ ~0.8fm. (38) 

The last result is an estimate based on the p-meson/V- meson radius-ratio [49, 74]. 

The Faddeev equation's kernel is completed by specifying that the quarks are dressed, with two 
of the three dressed-quarks correlated always as a colour-3 diquark. Binding is then effected by the 
iterated exchange of roles between the bystander and diquark-participant quarks. A Ward-Takahashi- 
identity-preserving electromagnetic current for the baryon thus constituted is subsequently derived [75]. 
It depends on the electromagnetic properties of the axial-vector diquark correlation. 



5.2 Nucleon form factors 

A study of the nucleon's mass and the effect on this of a pseudoscalar meson cloud are detailed in 
[76]. Lessons learnt were employed in a series of studies of nucleon properties, including form factors 
[77, 78, 79]. The calculated ratio fx p G p E (Q 2 ) /G P M (Q 2 ) passes through zero at Q 2 « 6.5 GeV 2 [78]. For 

2 

the neutron, in the neighbourhood of Q 2 = 0, fi n G r E {Q 2 ) / G r l i {Q 2 ) = —^Q 2 , where r n is the neutron's 
electric radius [79]. The evolution of /i p G p E (Q 2 )/G p M (Q 2 ) and H n G E (Q 2 ) / G n M (Q 2 ) on Q 2 > 2 GeV 2 are 
both primarily determined by the quark-core of the nucleon. While the proton ratio decreases uniformly 
on this domain [77, 78], the neutron ratio increases steadily until Q 2 ~ 8 GeV 2 [79]. 

A form factor ratio motivated by ideas from perturbative QCD [80] is depicted in Fig. 2. The 
parameter A is interpreted as a mass-scale that defines the upper-bound on the domain of so-called soft 
momenta in the perturbative analysis. A plausible value for such a quantity is A ~ Mn [77]. 14 It is 
interesting that the model of Ref. [77] yields neutron and proton ratios which cross at Q 2 ~ 8.5 GeV 2 . 
The model's prediction for truly asymptotic momenta is currently being explored. 

It is plausible that Nature's fundamental "constants" might actually exhibit spatial and temporal 
variation [81]. This motivates an exploration of the current-quark- mass-dependence of the nucleon 
magnetic moments, which complements work on hadron <r-terms [40, 82, 83]. Preliminary results for 
the quark-core contribution to this variation can be expressed through the following ratios evaluated at 
the physical current-quark mass [35]: 

— ^ = 0.016 — , — — = 0.0042 — . 39 

Hp m fi n m 

It is likely that pseudoscalar meson contributions will increase these values by a factor of > 10 [84]. 

The framework can naturally be applied to calculate weak and strong form factors of the nucleon. 
Preliminary studies of this type are reported in Refs. [85, 86, 87]. Such form factors are sensitive to 

14 NB. A value of A <~ 0.3 GeV corresponds to a length-scale r\ ~ 1 fm. It is not credible that perturbative QCD is 
applicable at ranges greater than the proton's radius. A reasonable value is A ~ 1 GeV, leading to riGcV ~ 0.2 fm. 
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Figure 2: Weighted nucleon Pauli/Dirac form factor ratio, calculated using the framework 
of Ref. [77] and presented with A = 0.94GeV: solid curve - neutron; dashed curve - proton. 
(Figure adapted from Ref. [35].) 



different aspects of quark-nuclear physics and should prove useful, e.g., in constraining coupled-channel 
models for medium-energy production reactions on the nucleon. This is important to the search for the 
so-called missing nucleon resonances 15 and the related problem of identifying exotic baryons. 

5.3 Quark orbital angular momentum 

Of significant interest is the distribution of an hadron's spin over the quark constituents and their 
angular momentum. In a Poincare covariant approach that can be calculated in any frame. The rest 
frame is physically most natural. The pion was considered in Ref. [39] and although J = the dressed- 
quarks carry significant orbital angular momentum. A mixing angle of ~ 30° can be attributed in the 
rest frame to the L = 1 components of the pion's Poincare covariant Bethe-Salpeter wave function. 

The answer is more complicated for the spin-| nucleon. In the truncation described above, a 
nucleon's Faddeev wave-function is expressed through eight scalar functions: no more are needed, no 
number fewer is complete. Two are associated with the + diquark correlation: Sip, and six with the 
1 + correlation: Ai,,,,^. In the rest frame in this basis one can derive the following "good" angular 
momentum and spin assignments, which add vectorially to give a J — | nucleon: 16 



L = 0,S= \ 


L = l,S=± 


L = l,S=l 


L = 2,S = I 


Sx , A 2 , Bi 


<S 2 , Ai , B 2 j 


c 2 





B x = IA 3 + IA 5 , B 2 = |A 4 + |A 6 , d = A 3 -A 5 , C 2 = A 4 -A 6 . (41) 

These assignments are straightforward to demonstrate and understand; e.g, in the rest frame of a 
relativistic constituent quark model the 5^2 terms correspond, respectively, to the upper and lower 
components of the nucleon's spinor. 

15 This refers to the fact that numerous excited states of the nucleon (N*), which are predicted by the ST/ (6) £g> 0(3) 
constituent quark model {CQMqxs), are not seen in the baryon spectrum that is obtained from amplitude analyses of 
nN elastic scattering. At this time more than half of the low-lying states predicted by CQMq X 3 are missing [88]. 

16 Equation (3.35) or Ref. [89] contradicts Fig. 6 of that reference. Equation (41) herein describes the correct assignments. 
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To exhibit the importance of the various L-S correlations within the nucleon's Faddeev wave- 
function, Ref. [62] reported the breakdown of contributions to the nucleon's canonical normalisation: 17 
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To illustrate how Eq. (42) should be understood, note that the largest single entry is associated with 
<Si <E> Si, which represents the quark outside the scalar diquark correlation carrying all the nucleon's 
spin. That is the w-quark in the proton. However, it is noteworthy that a contribution of similar 
magnitude is associated with the axial-vector diquark correlations, which express mixing between p- 
and ci- waves; viz., C\ <E> C 2 + C 2 <E> C\ . With C 2 all quark spins are aligned with that of the nucleon and the 
unit of angular momentum is opposed: (q |)(cjf 1)(P |), while with C\ all quark spins are opposed 
and the two units of angular momentum are aligned: (q l)(q l)(q l)(D j). This contribution is more 
important than those associated with S 2 ; namely, scalar diquark terms with the bystander quark's spin 
antiparallel. Finally, for the present, in this context one single number is perhaps most telling: the 
contribution to the normalisation from (L = 0) <8> (L = 0) terms is only 37% of the total. 



6 Coda 

The basic problem of hadron physics is to solve QCD. In order to achieve this goal a joint effort 
from experiment and theory is required. Our community now has a range of major facilities that are 
accumulating data, of unprecedented accuracy and precision, which pose important challenges for theory. 
It is the feedback between experiment and theory that leads most rapidly to progress in understanding. 
The opportunities for hadron physics will grow because upgraded and new facilities will appear on a 
five-to-ten-year time-scale. 

This volume provides an excellent illustration of the splendid experimental accomplishments of 
recent times, the achievements of theory in explaining and interpreting the observations, and some of 
the urgent open problems. Two of the things that stand out are the diversity of theoretical approaches 
and the interplay between them that is necessary to draw valid conclusions from contemporary data. 

One need not think too hard in order to list a number of the key open problems. It is essential to 
understand the origin and the nature of confinement. Asymptotic coloured states have not hitherto 
been observed, but is it a cardinal fact that they cannot? As described herein, confinement is related 
to the analytic structure of elementary n-point functions in QCD and especially to their properties at 
infrared momenta. Is it in the essence of non-Abelian gauge theories that they possess an Hilbert space 
of physical states that is reconstructed from only a small subset of all the theory's Schwinger functions? 
In progressing toward an answer to the question of confinement it will be necessary to map out the 
infrared behaviour of QCD's /3-function. The upgraded and future facilities will provide data which 
will guide that process. However, to make full use of that data, it will be necessary to have reliable 
Poincare covariant theoretical tools which enable the study of hadrons in the mass range 1-2 GeV. On 
this domain confinement and dynamical chiral symmetry breaking are both germane. 

17 The entry in location Si £g> Si indicates the integrated contribution associated with Sf. The entries are reweightcd 
such that the sum of the squares of the entries equals one. Positions without an entry are zero to two decimal places. 
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It is known that dynamical chiral symmetry breaking (DCSB); namely, the generation of mass from 
nothing, does take place in QCD. However, the origin of the interaction strength at infrared momenta 
that guarantees DCSB through the gap equation is unknown. This ties confinement to DCSB. Further 
to this connection, at nonzero temperature and chemical potential are the deconfinement and chiral 
symmetry restoration transitions simultaneous? If so, why, under what conditions, and also, perhaps, 
in which class of theories? 

It is important to understand the relationship between parton properties on the light-front and the 
rest frame structure of hadrons. This is a problem because, e.g., DCSB, an established keystone of 
low-energy QCD, has not been realised in the light-front formulation. The obstacle is the constraint 
k + := k° + k 3 > for massive quanta on the light front [90]. It is therefore impossible to make zero 
momentum Fock states that contain particles and hence the vacuum is trivial. Only the zero modes of 
light-front quantisation can dress the ground state but little progress has been made with understanding 
just how that might occur. Furthermore, parton distribution functions must be calculated in order to 
comprehend their content. Parametrisation is insufficient. It would be very interesting to know, e.g., 
how, if at all, the distribution functions of a Goldstone mode differ from those of other hadrons. 
Calculation [91] and experiment [92] together are needed to address that issue. 

A last question in this brief survey, but certainly not the least, how do the properties of hadrons 
change within nuclei, and within very dense systems? Absent an answer to this, one is not addressing the 
subject of nuclear physics. The experimental archetype for this problem is the EMC effect. Experiments 
have shown that in-nucleus parton distribution functions cannot simply be obtained by adding together 
the distributions within the constituent nucleons. A satisfactory understanding of this effect is still 
lacking, almost twenty-five years after its discovery. Within this theme one also asks for information on: 
the equation of state for dense matter; the properties of the hadron to quark-gluon phase transition; 
and the evolution between and nature of the states of matter within compact astrophysical objects. 
Plainly, in addressing this question one is embarking on a far reaching quest. 

There is ample evidence of the constructive feedback between theoretical approaches. For example, 
Dyson-Schwinger equation studies and effective field theory inform calculations within lattice-regularised 
QCD, and this interplay will become increasingly powerful as the parameters of lattice-QCD simulations 
come closer to the physical domain. Notwithstanding the investment in and progress that has been 
made with lattice-regularised QCD, for example, there is a pressing present and continuing need for 
model building and application in hadron physics. No single so-called ab-initio approach is currently 
applicable to all the phenomena attributed to QCD at all the length-scales that can be explored. 
Therefore, understanding and synthesising the wealth of extant and forthcoming data requires the 
continued development of models. Such tools answer an immediate need and provide a ready means 
to develop intuition about a complex system. Naturally, a reasonable model should obey all relevant 
constraints that ab-initio approaches and symmetries apply. That being so, it can then provide a 
rational and flexible guide in connection with our rapidly changing experimental environment. 

Continued investment in each of the nonperturbative tools that might provide insight into QCD 
is also necessary. Perturbative QCD is understood but its application begs a key question: at just 
which scale does its employment become valid? The answer can only be attained nonperturbatively. 
A number of experiments are described as exhibiting precocious scaling; i.e., they are associated with 
cross-sections that appear to evolve with momentum according to perturbative QCD expectations. 
However, only a comprehensive nonperturbative analysis can explain whether that is the result of an 
accidental cancellation over the - often small - momentum domain accessible to the experiment or a 
clear signal of perturbative QCD. Partial answers to this and related questions are possible through 
internally consistent model calculations that can unify many observables, e.g., Ref. [93]. But a complete 
answer will require truly nonperturbative calculations of the scale-dependence of observables within a 
framework that is able to exhibit the perturbative QCD limit, e.g., Ref. [39]. 

This contribution has outlined the utilisation of DSEs in studying some of the many facets of 
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the basic problem of hadron physics. The DSEs provide a natural framework within which light- 
quark confinement can be defined, and the origin and consequences of DCSB explored. To build 
an understanding of confinement, it is essential to work toward an accurate map of the interaction 
between light-quarks. Among the rewards are a precise connection between confinement and DCSB, an 
accounting of the distribution of mass within hadrons, and a realistic picture of hybrids and exotics. It 
cannot be too strongly emphasised that DCSB is remarkably effective at generating mass. For light- 
quarks it is far more important than the Higgs mechanism. It is understood via QCD's gap equation, 
which delivers a quark mass function with a momentum-dependence that connects the perturbative 
domain with the nonperturbative, constituent-quark domain. The existence of a sensible truncation 
scheme enables the proof of exact results using the DSEs. The scheme is also tractable, and hence the 
results can be illustrated and predictions made for observables. The consequent opportunities for rapid 
feedback between experiment and theory provide openings through which an intuitive understanding 
of numerous nonperturbative strong interaction phenomena can be reached via the DSEs. 
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